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FDulANhi^T  Shl-UFi.  fOti  Ib.J  CALCULATORS 


Introduction 

^  The  piiroose  of  thi-s^paoer  is  to  discuss  th®  efficiency  of  different 
methods  of  computation  wherein  the  elementr.ry  functions  and  othorr  must  be 
evaluated  on  the  IB’/  Fleet ronic  Calculating  Punch,  models  604  and  6C5,  and 
the  IBM  Card-Programmed  Electronic  Crdculator,  models  I  ai'.d  TI.  Thes^  mtrehMves 
widi-bo  referred  to  hereafter  simply  as -the  604  ,  605,  CPC-I,  and  CPC-ll, 
respectively.  Come  of  the  following  is  applicable  to  comnutation  cn  any 
electronic  computer,  but  those  computerr  writh  a  large  memory  cappcitv  and 
high-speed  internal  programmin£»  largely  remove  the  ver^  restrictions  which  make 
efficient  evaluation  of  functions  a  oroblem  on  the  604  cr  CPC. 

The  fact  that,  given  enough  time,  patience,  and  cards,  practically 
any  computation  can  be  performed  on  the  604  or  CPC  is  trivial.  What  is  of 
interest  is  how  any  job  from  a  fairly  large  class  can  be  most  efficiently  done 
in  the  day-by-day  work  of  cn  IP’l  comouting  establishment.  Assu.minr  that  a  prob¬ 
lem  is  already  stated  explicitly  -  i.e.,  that  the  computing  group  is  required  to 
evaluate  certain  items  or  obtain  numerical  solutions  to  certain  equations,  or 
both,  for  a  given  numt'cr  cf  parameter  combinations  -  it  would  still  be  necessary 
to  make  several  decisions  and  to  perform  considerable  preparat ory  work  before 
putting  the  job  on  a  machine.  Hence  it  is  aonarent  that  standarc'i  z®d  methods  Cenn 
be  as  much  of  an  aid  to  computation  as  the  actual  ■machines,  provided  t^  a t  such 
methods  are  general  enough  to  embrace  a  vrlde  variety  cf  Jobs  and  ar®  not  cumber¬ 
some.  Also  It  is  helpful  to  have  set-ups  available  which  allow  decisions  to  be 
made  on  the  basis  of  the  problem  presented  without  excessive  consideration  of 
how  such  and  such  a  step  of  the  problem  can  actually  bo  calculated  efficiently* 
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Floating?  Decinal  (FD)  S(^-l)ps  for  the  CPC 

For  the  CPC,  gener.'il-purpor.e  set-uos  (or  coeci al-nurpo^e  s^t-uns  for 
a  lon^  continuing  orohlrin)  are  practicallv  a  necessity  if  the  ti’ne  reculrcr* 
for  coding,  wiring,  rnri  testing  is  to  be  kept  within  reisonalle  limits.  The 
CPC  may  be  looked  upon,  not  as  a  computer,  but  as  r  neat  oackare  of  oarts 
which  the  set-up  man  assembles,  bv  means  of  vrlres,  into  a  coT.outing  system. 

The  skill  which  the  coders  then  acquire  in  using  this  system  becomes  one  of 
the  assets  of  the  organization.  If  the  ryctein  is  complicated  to  use,  fewer 
coders  will  use  it  or  become  proficient  with  it.  If  the  system  is  elaborate 
enough  to  be  truly  general  and  still  simple  to  use,  it  is  almor*-  bound  to  be 
complicated  plug-board-wise;  hut  if  the  speed  of  oneration  is  thereby  reduced, 
the  advantage  of  simnlicity  in  preparing  a  job  is  If-rgely  overcome  by  the  t^me 
consumed  in  running  it.  This  situation  is  aggraivated  by  flo'-ting  decimal  (FT) 
operation,  which  many  peoule  pr-efer  for  the  CPC  rand  which  certainly  has  advantsf^es 
if  it  is  properly  used.  The  writer  has  spent  conslderahle  time  in  im-'rovlnm 
and  developing  general-purpose  FD  set-ups,  first  fer  the  CPC-I  and  then  for  the 
CPC-II.  Mnny  of  the  techniques  used  are  aim,  r.oDli  cable  to  fixed-no^nt  rn'rallel 
work  on  the  604. 

A  general-purpose  FL  set-uo  for  tf'e  CPC-I  can  include  square  root 
along  with  the  four  basic  aritlimetic  anc  necessary  rr.Jtninulr t1  ve  ooer-atiens. 
Pingle-card  evaluation  of  other  functions  ’s  virtually  orecluded  for  8xfl  arith¬ 
metic  by  the  limitations  of  the  mrichlrie.  To  fill  this  need,  standardized  sub¬ 
routines  are  often  orovlaed,  many  using  rational  anuroxl'’ationF  to  the  functions. 
Mr.  Cecil  Hastings,  Jr.  has  develonrd  many  vsiluable  arnroxi  matl  nns  of  this  kind.^^ 
Such  approximations  are  equally  valuable  in  604  work,  where  they  ar»  incorporated 
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in  the  procedure  r&ther  than  exirtinr  hs  sub-routines.  An  inherent  riavnntnpc 
of  these  sub-routines  Is  that  the  lenfrth  of  the  comoutation  can  be  readilv 
fitted  to  the  accuracy  desired  by  insuectini?  the  sfiecificati or.s  of  the  auoroxl- 
raation,  whereas  the  number  of  terns  required  in  f  series,  for  examnie,  will 
vary  widely  for  a  moderate  interval  for  the  rmiment. 

The  CPC-II  makes  it  oossible  to  consider  FD  set-uos  of  a  much  more 
elaborate  kind.  Its  hi/^her  cost  1s  worth  v.hile  cnlv  If  greater  efficiency  can 
be  obtained.  This  is  not  as  easy  to  do  for  a  general-puroose  set-up  ar  first 
appears,  since  the  internal  storage  of  the  605  is  no  gre<ter  thaji  th;  t  of  the 
604  and  there  is  actually  less  "compute  tine"  available  ncr  card  cvcle  on  the 
605.  At  least  three  v.ays  to  increase  the  output  of  useful  work  are  as  follows: 

(1)  to  perform  more  ooeratlons  oer  card  in  the  same  time, 

(r)  to  include  functions  ns  single-card  ocerftionr  instead  of 

u  si  ncr  sub-rouM  nes,*  and 
(3)  to  carry  more  digits  for  each  number. 

The  last  item  mentioned  is  a  subject  of  its  own  and  will  not  be  of  concern  here 
except  in  what  immediatelv  follows.  The  size  oi'  the  storage  units  v.ithin  the 
604  or  605  limits  the  programmer  to  8x8  FI-  arithmetic  if  any  roced  is  tc  be 
maintained  end  if  very  many  operations  are  to  be  built  into  the  calculator  board. 
Furthermore,  since  the  main  storage  for  the  CPC  is  in  the  941*  t  ,  aval  la’' le  only 
via  10-digit  channels,  and  since  it,  is  necessary  to  carry  a.n  exponent  of  10  for 
tiie  scale  factor  (usually,  although  not  vpw  accurately,  called  a  "power")  with 
each  number,  8x8  arithmetic  with  7-digit  pov.err  seems  to  fit  the  CPC  t-et.ter  than 
♦Increased  sneed  is  of  course  sought,  but  even  without  it,  the  elimination  of  sub¬ 
routines  is  an  advantage  by  saving  storage  units  in  the  tab  or  d/,1 ,  cooinr  and  key¬ 
punching  or  else  reproducing  time,  and  card  handling. 
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any  other  combination  for  peneral-purpo^e  FD  cet-ups.  To  carry  more  dirits, 
even  10  or  12,  some  severe  restrictions  must  be  olaced  on  Hcidrersei  -^f  the  oper¬ 
ands  and  allowable  operations.  Double  precision  boards  (13  or  ’^ore  dipits  nor 
number)  are  of  course  voiy  useful  for  some  r.ooli  cat  ions,  and  the  CPC-IT  nrkor. 
it  possible  to  have  much  better  ones  than  the  CPC-T  allows,  hut  they  prohibit 
sinple-card  evaluation  of  functions,  excentinr  nerhrrs  scusre  root. 

A  Four-Address  FD  Set-Up 

On  the  other  hand,  the  only  ncssible  sdvanta<je  of  FT  ooeranrir  v.ith 
fewer  than  8  dirfts  is  to  reduce  conputin^^  time,  since  the  channels  and  the  041 
units  cannot  be  adapted  to  efficient  hanclinr  of  "^^-words."  Hut  the  ^-05  v.ill 
go  through  60  program  stens  every  card  cycle  anywpv,  sc  that  if  2  ariitrary  FD 
operations*  on  8-digit  operands  can  be  nerforned  In  these  60  steps**  (without  an 
excessive  number  of  multi  plications  or  divisions)  there  i s  no  reason  to  reduce 
the  word  length.  Iv.  is  couarent  that  addition  is  the  rent  difficult  arithmetic 
operation  '.vith  FD  numbers  since  it  involves,  rctuallv  four  operations: 

(1)  a.  right  shift  to  align  the  decimal  noints  of  tho  adnends, 

(2)  the  addition, 

(3)  the  checking  of  the  sum  for  a  carr/  or  cancel  lati  on  of  lescinp 
digits,  with  approoriate  right  or  left  shifts,  and 

(4)  the  correction  of  the  cower  cf  the  sum. 

Although  these  different  parts  of  the  addition  routine  may  Vjp  used  an 
adjuncts  to  other  opf^rations,  th^^  complete  routine,  nevert.he]  r  s,  consumes  over 

*More  than  2  such  operations,  undulv  comnlicrte  not  or.lv  vnriiig  but.  also  coding. 
**A  single  "program  repeat"  will  always  cause  a  delay  cycle,  even  without  using 
"repeat  delay";  at  least  this  is  true  on  tiie  three  CPC-II's  which  the  vaiter 
has  used. 
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half  a  sweep  -  43  progr&n  steps  being  the  best  the  writer  has  ijeen  3*^10  to  do 
to  allow  for  all  eventualities,  although  2  or  3  of  these  steps  may  have  been 
saved  In  combination  operations.  Thus  If  2  arbitrary  operations  are  to  be  per¬ 
formed  in  one  sweep,  two  banks  (chains)  of  programs  on  the  605  must  be  used  Just 
for  the  3  basic  arithmetic  ooeratlons  ( subtrcicti  on  can  be  accomplished  by  changing 
the  sign  of  the  subtrahend  on  the  tab  board),  one  bank  performing  addition  first, 
followed  by  multiplication  or  division,  and  the  other  bank  doinm  the  addition  last. 
The  writer  has  wired  such  a  set-up,  which  takes  tiiree  8-di  r*!  t -pi us-2-di git -power 
inputs  (cny  two  such  inputs  from  the  941* s)  and  performs  any  possible  combination 
of  2  of  the  4  basic  arithmetic  operations  on  each  card  cvclc  at  150  cards  per 
minute,  except  for  2  additions  (or  subtractions).  The  latter  co'rhinrtion  is 
allowed,  but  a  delay  is  required.  (Occasionally,  other  operations  require  a 
delay,  when  there  are  many  8’s  or  9's  in  a  multiplier  or  quotient.)  This  set-up 
is  proving  t-  be  very  useful,  since  it  also  allows  computation  of  square  root, 

10*,  log,.x,  and  cos  x  as  single-card  operations.  HI  of  these  ooerati'^ns  are 
quite  fast  except  log  x,  which  is  suVi-prcgra.Timed  (Special  Progra-mminp''  on  the 
tabulator).  The  only  interval  restriction  on  arguments  is  f^r  cos  x,  for  which 
only  0  X  ^  TT  is  allowed.  Series  ere  used  for  10^  and  cos  x  ajid  an  rncroxi mrti  on( 
Is  used  for  log  x.  Since  rcort  computations  consist  of  more  basic  nritlunetic  op  r~ 
ations  than  anything  else,  and  hence  the  faster  they  can  be  done  better,  this 
set  of  functions  and  operations  provides  a  pood  compromise  for  increasin  ’  the 
useful  output  of  a  CPC-II  by  the  first  two  means  stated  previously. 

An  FD  Set-Up  With  Several  Functions 

For  jobs  requiring  frequent  use  of  funct^’ons,  cne  mav  well  be  satisfied 
with  single  arithmetic  operations  if  a  greater  vrrietv  of  functiens  can  be 
evaluated  at  electronic  speed  as  single-card  operations.  The  RAI^L  Computing 
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Group  has  one  such  set-up  (FD)  for  the  CPC-II 
card  operations  are  allowed,  all  of  which  are 
605: 

1.  A  +  B 

2.  A  -  B 

3.  AE 

4.  A/E 

5.  sin  E,  -  H’f.  X  <  Tr 

6.  CCS  B,  0  £  X  £:  TT 

7.  10^,  -51  i:  X  4  49 

9.  VjT,  10“^^  4X4  10^"^ 

11.  AB  +  P 


on  which  the  followln/'  single- 
evaluated  conpletely  within  the 

1. '  -A-B 

2. '  5  -  A 

3.  '  -AB 

4. '  -A/B 

5. '  sinh  E,  -2.9  <  x  <  2.Q,  approx. 

6. '  cosh  D,  04X4  2.Q,  nnnrcx. 

8.  Io^t^qX,  10-51  ^  <  ^0^^ 

10.  fircsin  x,  -1  £  x  <1  (|  x|>  1 
treated  cc  |  x  |  -  1) 

12.  AP  +  B 


A  and  B  are  arVitrary  operands  and  P  rterns  the  orevious  reruit,  although  as  far 
as  the  605  is  concerned  it  could  he  arbitrer'''^.  Series  nre  used  for  all  the 
functions  except  -y/lT  ,  for  which  the  usual  Newton's  Ticthod  is  used,  2  iterations 
being  performed  on  one  sweep.  Fpeed  and  scciiracy  are  both  fairly  rood  for  the 
functions.  Thirty  arguments  between  0.1  and  1.0  for  arcrin  x  were  run  in  one 

n 

minute  with  a  maximum  error  of  3  ^  10“  .  An  interesting  si  denote  is  that  it 
took  some  hunting  to  find  the  proper  argument  which  v.'ould  mive  this  maximum  error. 
It  depends  on  the  relative  size  of  the  rounding  error  in  takinc'  x^  to  8  decimals 
os  compared  to  the  size  of  Vl-^^  ,  since  for  oranments^  greeter  than  0.7,  the 
relationship  arcsin  x  =  ^/2  -  arcsin  VT-^  is  used.  (Only  non-negative  argum^'nts 
are  used  in  the  series,  th(f  sign  being  handled  seonratel  y. )  The  argun'^nt  0.^‘^'^O3 


gave  the  maximum  error. 
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FD  Techniques  In  Fixed-Decimal  Computation 

All  the  above  functions  exccot  lop  x  and  '/jP  rcGi'ire  tUt  the 

argument  be  put  in  essentially  fixed-decimal  forr.  rithin  the  605;  thus  they  can 

be  evaluated  just  as  well  on  a  604  for  parallel  v;ork  or  on  a  CPC-I  for  fixed- 

point  arithmetic.  However,  Fb  notation,  at  least  within  th»^  605,  is  essential 

for  log  X  and  a  decided  advantage  for  \Jir  .  To  evaluate  them  efficiently  on  a 

604,  a  ^'time  emitter  is  needed  on  the  521  sc  that  th^*  argument  cm  be  left- 

shifted  on  read-in  until  there  is  a  leading  non-zero  di  pit  for  tlK'  lor  oo^ration 

or  a  leading  pair  of  digits,,  not  both  zero,  for  the  souare  root  oD^ration,  the 

shifted  argument  then  being  considered  as  a  decimal  number. 

The  advantage  of  1  eft-shi ftinp  the  argument  before  anilymp  Newton's 

method  to  square  r*ot  Is  cue  to  the  fact  that  if  the  fii-rt  m-rrs  is  much  larger 

than  the  true  root,  many  iterations  are  requirec  before  the  first  (non-zero) 

significant  figure  is  obtained.  The  most  convenient  starting  value  f'-'p  the 

whole  interval  .01  <  x  <  1.0  (which  allows -for  thf’  fact  that  x  can  onlv  be 

shifted  an  even  n\jmb<^r  v-^f  places)  is  Tn  fact,  for  x  _>  .onnqn^ 

3^x  +  =  \/x  to  8  places,  ^  fact  useful  in  handling  the  troublesome  st  <^cirl 

case  of  X  =  .99999999,  which  causes  an  overflow  on  =  ^(Vq 

other  cases,  y^  will  be  larger  than  either  x  or  /IT  ,  fiid  convergence  will 

proceed  with  no  complications.  Of  course,  the  root  thus  found  mu^t  be  rirbt- 

rhifted  back  half  the  number  of  olrce*-'  the  ar>mjmpnt  ypt  1  eft-sbi  ft ed.  Tn  FP 

set-ups  the  power  merely  h' s  to  be  corrected,  an  ooeraticn  ersilv  serformed 

by  dividing  the  original  Dov.»^r  olus  50  by  two  end  adding  the  reT;e’'drr  to  the 

quotient.*  This  is  done  first  and  if  t.hc  r'^mainder  is  non-z-^ro  (i.e.,  l). 

found 

then  the  argument  itself  is  ri  giit-shi  fted  one  olace.  The  root/vill  alwavs  lie 


*This  assumes  unity  is  represented  by  10000000  51. 


p-m 

prif'^o 

in  FP  form  iP  thr*  r.rgrinrnt  was  orlpinallv  fnri  will  bf‘  corroct  In  pv.y  cnr.c. 

Series  f'>r  the  Uementary  Functions 

It  i  s  convenient  to  divide  the  power  series  for  the  el^'nent'rv 

1  unctions,  for  purposes  of  m-chine  coifiputsticn,  into  three  cles.sts.  The 
00 

exppjision  e^  =  2:  x"/n!  converger  quicklv  to  p  soec^’ fi  ed  .'■•ccurr.cy  because  of 
o 

the  n!  in  the  denoTiinntor  of  the  tern  t^^.  Thus  the  ratio  r(n)  =  i^n+1  /  “n  h^ 

x|yn+l,  r.nc  as  soon  as  n  >  k  x  ,  r(n)  <.  l/k  and  the  succeedin;^  t‘=rns  quickly 

decrease  in  magnitude  if,  say,  k  >  2.  For  small  intervals,  sey  |x  c  iof»^10<  2.t, 

n  =  5>  2  X  ,  Also  t  <  (.'.5)V5!  <  1  snd  t  !<  2“^.  Anv  series  of  the  form 

5  5+m 

y~  whore  c^  is  of  the  order  of  l/n!,  will  he  called,  nriitrarilv,  exponent  ^  al- 

type  series,  for  want  of  a  Irttcr  name.  They  converge  uniformly  over  ru-iy  interval. 

They  include  e^,  sin  x,  cos  x,  sinh  x,  and  cosh  x. 

If  c^  in  the  above  series  is  of  the  order  of  l/n,  thr  *  erics  will  be 

called  an  harmoni c  type.  This,  tyj^e  induces  the  serie?  ^or  log  x,  tan“^x,  and 

tantr^x. 

ftill  a  tided  type  of  Eeriei  is  met  with  in  the  oiRT,rntar,r  Air.ctions 
but,  for  practical  purposes,  may  i)e  cnns1der»=d  alrnm  v.iti  tt  o  *aro:onic  tvoe.  Th^  s 
t^Toe  consists  of  the  seri'-s  who 
The  first  fraction  may  be  wr'tt 

u-'l  /T\ 

00  ^  00 

all  n.  In  fact,  since  ^  1_  diverges^,  t  }v  infinite  nroc.uct  TT  (1-  1_)  divrrr'*s 

1  2k  2k' 

also  (in  the  stricter  sense)  and  h^nce  diverge;  to  zero,  decrcrsUiP  monot callv 
a;;  n  increases.  Therefore  these  orociuct  coefficientr  cann''t  hurt  tii'’  cerverrence 
of  the  series  in  which  they  appear  nor  cm  t’iPy  inerrr  so  th*^  si  so  of  t^o  neglected 
terms  when  tiie  series  are  tmncrt.ed.  On  tl”  ot^'m  hand,  they  do  ni't  help  c-^nvergence 


se  n^^  term  i  s  of  the  form  1  •  ^ ■  5»  • « ( 2n-l ) » 

n  r*/!.*6»*"(2n)  2n+l 

I  f  (T-  which  is  obviously  less  than  1  ^ci 


very  much,  at  least  I'cr  intervals  oracti  caide  for  cc^outati  on.  It  is  tnie  that 
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+  V  l*3«5‘-(?n-l)  .  X 
Z_  2*4*6*  ••(2n)  2 


•*•  =  a  resin  x  convorercs  for  x  =  1,  whereas 

2n+l 


00 

r 


^n+1  1 

_ =  tanh“-‘-x  does  not,  but  the  convergence  of  the  first  series  is  too 

2n+l 


alow  for  X  near  1.  It  converges  to  a  specified  accuracy  in  a  reasonable 
number  of  terns  for,  sny,  j  x  <0.7,  l^Jt  to  does  the  second  series.  The  only 
other  series  similar  to  that  for  arcsine  occurring  among  the  elementary  functions 
is  the  series  for  sinh“lx,  which  is  the  alternating  series  of  the  same  form.  Such 
series  with  even,  insteao  of  odd,  tonr.s  are  esnentiadly  polynomials  in  arcsin  x  or 
sinh""^x.  However,  some  other  functions  can  expressed  in  similar  forms  and 
evaluated  in  the  604,  e.g. ,  complete  elliptic  integrals  of  the  secono  kind, 

E(k)  =  J  (l-k2Eln20)idei.  Complete  elliptic  integrals  of  the  first  kind  cm  he 
easily  expressed  as  double  infinite  series,  nf  a  similar  nature,  but  the  writer,  at 
least,  has  not  yet  been  able  to  adapt  them  to  604  computation  as  direct  operationt 
because  of  the  difficulty  in  summing  the  coefficients  of  each  nower  of  k  or, 
alternatively,  summing  the  powers  for  eacii  coefficient.  The  series  for  F(k)  is 
discussed  subsequently. 

Evaluating  The  Log  Function 

The  most  useful,  harraonl c-tyoe  series  are  those  for  log  x.  Here  a 
leading  non-zero  digit  is  not  merely  e  speed-u”  device  but  essential  for 
practical  convergence  and,  in  fact,  to  the  use  of  eporoxi mat  i  ons  as  well,  sir  ce 


the  latter  are  usually  del'ined  for 


Ql.l-o] 


or  some  interval  wh^‘ ch  may  easily 


be  transformed  into  this  one.  For  fixed-ooint  arithmetic,  the  number  of  places 
the  argument  is  left-shifted  determines  the  characteristic  of  the  common  log;  in 
FD  work  the  f)Ower  amounts  to  the  seme  thing.  The  natural-log  series  is  easily 
adapted  to  taking  common  logs  by  using  log^QX  -  M’log^x,  where  M  =  -  .4'’/.?'^44fi. 

The  advantage  of  taking  common  logs  is  of  course  in  the  use  cf  the  characteristic, 
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since  this  is  independent  oT  the  EiCturJ  dibits  in  the  nunher,  clfce  zero's  riis- 
reperded. 

00 

The  series  5(>)  =  ^  jP/n  =  lop  (1+x)  -  lop  t  is  convenient 

1 

for  eve.luttinp  lops.  This  series  only  ccnverres  for  -1  <  x  <1,  but  ever;  this 

interval  is  too  1;  rpe  to  ie  of  urfi-ctlcel  interest  in  corrut st  < cn.  For  x  =  1,  n 

hundred  million  ter^.s  woulc  be  needed  to  keep  the  error  lees  tbrn  IC'^I  For 

-14  X  <  0,  S(x)  is  not  even  rltenietinp,  and  even  though  the  reries  converpes 

»tnd  r(n)  <.  1  for  all  n,  it  cannot  be  puaranteed  that  the  errnjr  cerr-Titted  by 

truncating  the  series  is  nc  prerter  thfn  the  -napnltude  of  the  first  term  ''matted. 

Fortunately,  hoaeverl  for  x  4^  .  S  Kn)  ^  .5;  in  this  esse  the  error  comitted 

by  truncatinp  the  series  is  less  than  t  ^2”  =  t  ,  where  t  is  the 

I  m  I  I  ra  I  m 

last  term  includec,  and  this  is  sufficient  for  the  onniose  at  hand  If  m  <8  not 

too  large.  If  |  ^„.  |  4.  10"®  is  required,  it  would  mean  tift,  for  jxj  /  .5,  .5®/m<10~®. 

A  simole  calculation  shows  that  this  is  satisfied  for  t.  =  21,  which  in  a  quite 

modest  number  of  terras  lor  electronic  computation.  Restr'etinp  x  to  .5) 

does  not  render  the  series  useless  by  any  means,  in  fact  it  is  quite  fdeoucte  for 

fvaluatinp  lop  x  for  a  very  larpe  interval  if  the  FP  notation  if  used  and  if 

proper  but  very  simple  transforr.atlcns  are  made.  In  contra.'t,  the  interval  f-.S  .5) 

is  oracticelly  useless  for  the  exponent  i  al-tvoe  series.  Thus  the  d^  r,advoijtape  of 

the  slow  convergence  of  S(x)  is  larpely  overcome  bv  the  small  interval  needed. 

Keeping  x  in  (-.5,. 5)  requires  t  in  (.5,1.5).  Vultiplvinp  the  first 

term  of  the  series,  t, ,  by  U  and  develooinp  t  from  t  transfoms  the  natural 

X  n^i  n 

lop  to  the  coomicn  lop,  and  ti:i  s  can  onlv  h^ip  converpence.  The  FD  notation,  or 

the  equivalent  left-shi  ftinp,  automatically  restricts  t  (i  priori  np  th«  characteristic) 

to  [~.1,1.0)  unless  t  =  0;  special  nrovisions  can  be  made  to  detect  inadvertent 

zero  arguments.  To  raise  the  lower  bound  on  t,  the  relationship  lor  t  -  lop  kt  - 

log  k  can  be  used.  It  is  convenient  to  use  a  power  of  2  for  k,  since  “  .'’CIOXXX) 
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and  hence,  even  if  this  number  must  be  emitted  during  electronic  programming, 
it  does  not  take  many  programs  because  there  are  only  3  non-zero  digits.  On  the 
605,  both  3’ a  can  be  emitted  at  once.  By  testlnp  the  lending  dipit  of  t  with  a 
digit  selector  at  read-ln,  kt  con  be  made  to  lie  in  Q6,1.6)  as  foliowo: 

Leading  digit  of  t  k  kt  in; 


0 

1 

2 

3»4,  or  5 
7,8,  or  9 


?JtR0R. 


8 

4 

2 

1 


[^8,1.6) 

[8,1.2) 

.6,1.2) 

[.7,1.0) 


[• 


If  t  is  slightly  less  than  .2,  kt  =  1.6-;  but  then  x  =  kt-1  =  .6-  nnd  the  series  is 

alternating,  which  improves  the  accuracy,  for  now  the  truncation  error  is  no 

,  n  ^  -8 

greater  than  the  first  term  omitted.  For  x  =  .6,  if  t^^  =  (.6  /n)  10  is 

required,  n  must  he  allowed  to  reach  29  -  certainly  not  too  many  terms  to  take 
for  the  extreme  case. 

In  the  second  FD  set-up  described  ah^ve,  nn  even  simpler  choice  v.-an 
made  for  k,  namely,  k  =  4  if  kt  4  1.4,  k  =  1  otherwise.  This  only  ruarantees  kt 
to  lie  in  [.35,1.4),  but  verv  good  results  are  obtained  nevertheless.  Since 
logi(^t  i  E  much  less  than  1  in  this  interval,  8  decimals  are  crrried  throughout 
the  series,  the  correction  -. 60206  =  log  I/4  being  added,  if  necessary,  at  the  end 
^to  prevent  overflow).  Finally,  the  oroper  characteristic  is  added  and  the 
results  rounded  to  7  decimals,  which  are  usually  good  all  the  wav  out. 

A  fixed-point  604  board  was  wired,  using  the  left-shifting  technique 
and  the  above  table  for  k  and  evrJ.uat1np  2  terms  per  sweeo.  IJsinp  rale  selectors 
to  move  the  (actual)  decimal  point,  logs  could  be  evaluated  to  7  good  dec  Vais  for 
arguments  lying  between  10"^  iuid  1C? -10“^.  However,  the  speed  wts  between 
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33^  and  50  cards  per  minute.  A  faster  board,  practicollv  a  constant  50  cards 

per  minute,  is  now  used  which  pives  virtually  the  same  accuracy.  This  604  set-up 

(2) 

uses  a  rational  aoproximation  of  the  form 

lop  t  =  f(x)  =  +  Cr^x^  +  C^x^  +  +  C^x  +  .5, 

where  the  C.  •  s  are  9  di^it  decimals  and  x  *  jL-Z — for  1  ^  t  £  10.  IJsinp 
^  t  +  yib“ 

X  =  f(x)  -1,  instead,  allows  t  to  lie  in  f. 1,1.61.  The  Ci '  c  are 

i(t  +  vCl)  •- 

rounded  to  8  clpits;  anci  t  (8  dipits),  v^K  to  8  cipits),  C^,  and  Cj  are  fed  in. 

C^,  n  =  5»  3»  1»  is  developed  vdthin  the  604  bv  multiplyinc'  ^^^2  emitt.inr  proner 

corrections.  This  board  was  much  more  difficult  to  propram  than  the  series  hoard, 
but  it  is  faster.  For  CPC  boards  the  nnproximation  is  not  feariVle  v.ithin  the  605, 
since  it  would  tie  up  too  many  oropraPiS  for  one  function;  but  it  is  easily  sub¬ 
programmed  as  in  the  first  CPC-II  set-u[)  described  previously. 

Evaluating  Functions  V.lth  Exponent! al-Tvpe  Series 

The  accuracy  of  the  exponent! al-tyoe  functif.ns  is  limited  for  peneral- 

puroose  8-digit  set-ups  by  the  fact  that  x  must  be  allowed  to  exceed  1  if  the 

evaluation  is  to  be  very  useful,  Althouph  this  does  not  'mnecir  conv^rp'’'nce  and 
only  moderately  increases  the  number  of  terms  required,  it  doer  allow  the  individual 
terms  to  become  greater  thp-n  1  even  thouph  the  sum  mav  be  quite  small.  Thus  fever 
significant  figures  can  be  retained  on  a  heavily  nrogrfin.med  beard. 

To  use  e^  to  evaluate  10^  in  a  m.-jiner  somewhat  the  reverse  to  that  used 
for  logpt  to  obta.in  log^Qt,  x  must  be  allowed  to  varv  over  (-r.X)?5851,  2.'50?5851), 
i.e.,  (logg.l,  logglO),  since  the  decimal  cart  of  t,  d(t),  will  vary  ever  (-1,1)  and 
hence  lOd(t) 

over  (.1,10).  The  whole-number  -'art  of  the  argument,  *  can  of  course 
be  handled  separately  to  shift  the  decimal  point.  For  ppn<^ral  term  of  the  series 
for  e^,  jx*’*/'^!  ^  10  for  x  C  3.91Z.7  and  for  all  n,  in  oarticular  n  =  3,  since  the 

maximum  terra  for  3  <  x  <4  is  t^.  Hence  it  will  be  sufficient  to  allow  onlv  one 
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.ace 


digit  for  for  all  n  v.ith  QxJ  <  log^lO.  Hov^ever,  since  the  last  nl) 
is  bound  to  be  Inaccumte,  more  than  6  good  decimals  cannot  be  oxoccted  for 
e?^  or  10^  if  only  8  places  are  carried.  Similar  considerations  hold  for  sin  x 
and  cos  x.  In  fact,  it  would  be  desirable  to  allow  x  to  vary  over  ^-PV'  ,  P7rJ 
for  sin  x  and  cos  x,  but  this  interval  is  too  large  even  when  allowinp"  a  place 


for 


In  the  set-ups  previously  described,  a  reduction  of  the  argument 


to[-TT,'TrJ  must  be  card-programmed.  For  a  604  board  for  only  one  or  two 
functions,  this  reduction  could  no  doubt  be  programmed  internally.  The  argument 
for  sinh  x  or  cosh  x  is  restricted  by  the  range  of  the  function  rather  than  by 
the  magnitude  of  the  terns.  Again,  for  special  604  boards,  one  could  arrange  to 
store  the  tens  and  higher-order  digits  of  the  sum  separately,  thus  extending  the 
Interval  up  to  about  (-'^.9,  3«9)  while  keeping  7  decimals.  These  extra  devices 
Impose  too  great  a  load  on  the  programming  of  a  general-purnose  FD  CPC  set-up. 

For  some  jobs,  the  relative  error  of  sin  x  or  cos  x  as  comoared  w^’th 
X  Is  the  important  thing.  In  such  eases,  series  as  used  above  -  i.e.,  with  7 
fixed  decimal  places  -  are  not  of  much  value.  Mr.  Hastings'  nDproxinati ons  for 
sin  X,  as  well  as  some  other  functions,  give  a  rela.tive  error  rather  than  a.n 
absolute  one.  The  computation  can  be  carried  out  as  an  FD  routine  or  on  a  double¬ 
precision  set-up,  and  good  accuracy  could  be  obtained  for  verv  small  arguments. 
This  is  an  invaluable  method  for  such  Jobs. 

The  Arctangent  Function 

Perhaps  the  most  difficult  elerrientary  function,  normally  used,  to 
handle  by  series  over  an  adequate  interval  is  tan”^x.  The  series  for  this 
miserable  function  has  all  the  disadvantages  of  the  harmonic-type  series  with 
none  of  the  saving  graces  of  the  log  function.  The  first  complication  is  due 


to  the  fact  that  T(x)  =  ^  (-1)*^  x^^^  =  tan"^ 


-1.  = 


X  converges  only  for  j  x  |  ^  1. 
-1/ 


For| X  I  ^  1,  the  relationship  tan"^x  =  iTTCsgn  x)  -  tan”^(l/x)  must  be  used^ 
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r  1  if  X  >  0 

where  spn  x  “  if  x<  0.  fuch  a  choice  is  difficult  enough  to  propram  in 

itself,  and  there  is  the  even  greater  difficulty  that  the  convergence  of  T(x)  for 
I  X  near  1  is  theoretical  only,  since  the  order  of  :nagnituce  of  n  to  ensure  an 
error  of  less  than  10  °  for  such  values  of  x  is  10^.  T(x)  can  be  evaluated  on  the 
605  fo  that  jt^j  becomes  less  thrin  10~®  for  2n+l  <  100  if  ^|  <  -^709.  By  letting 
r(n)  =  mx  ^/m+2,  ra  -  fn+l,  only  the  odd-power  ter-ms  need  he  evaluated,  i.e.,  actually 
less  than  50  terns.  An  earlier  FD  C  P  C-II  set-up  than  those  previously  described 
included  this  function  instead  of  arcsine.  Fvexy  irnapinable  kind  of  reduction 
formula  was  tried,  but  all  proved  to  be  too  involved  to  program  intemailv.  Resort¬ 
ing  to  sub-programming,  the  followinp  procedure  was  finallv  used; 

Taking  the  well-known  formula,  tan  20  =  ,  let  x  =  tan  6,  y  = 

1-tan^ 

I  ej  <  irr,  whence  xy^  +  2y  -  x  =  0,  giving  y  =  -l^>/i-x^\  Since  |  v|  4.|  x  I  ,  the 

X 

plus  sign  must  he  chosen.  Applying  this  apain,  t  =  1+  \/l-y^  ~  e/4.  Fince  t  5r 

y 

.41421  =  tan"^B,  the  series  T(t)  will  converge  quickly.  Finally,  tan“^x  =  4  tan~^t. 

This  sub-programmed  reduction  introduced  very  little  loss  in  accuracy,  but  it  was 

slow.  However,  it  would  hsjxdle  any  argument  up  to  about  ^  .3I6  K  IC^^.  A  rational 

(3) 

approximation  would  probably  do  as  well  with  about  the  same  speed. 

The  arctangent  function  is  one  for  which  a  continued  fraction  should 
prove  very  useful  for  special-purpose  boards.  Since  tan'^x  is  an  odd  function,  the 
sign  can  be  "remembered"  with  selectors  on  the  tab  or  521.  The  epd  x  in  the 
transformation  previously  mentioned  may  th^n  be  ignored,  and  all  arguments  reduced 
to  [0.1]  by  not  more  than  one  division.  The  continued  fraction  for  tan”^x, 

tan~lx  =  X _ 

1  +  1 

3  +  4  x^ 

5  +  9  x^ _ 


7  +  . 
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is  valid  for  all  real  x  but  would  be  unwieldy  for  larpe  x.  It  it  unlikely  that 
the  maf^itude  of  x  could  be  tested,  the  continued  fraction  evaluated,  end  the 
correction  of  j'H'  -  tan-l(l/x)  made  when  necessans  all  within  a  ^04,  Hit  it  mi  pht 
well  be  done  in  2  operations.  It  is  impossible  to  evaluate  a  continued  fraction 
from  the  top  down  with  8-diirtt  numbers  in  a  fOA  or  605  because  of  stcrav**  limi¬ 
tations.  In  fact  it  requires  all  the  ftornge  to  evaluate  then  from  the  bottom  up. 
This  puts  continued  fractions  at  a  disadvantage  with  rescect  to  t ran sfo mat  ions  as 
compared  with  series.  The  continued  fraction  has  no  sum  to  which  a  correction  can 
be  pre-added,  and  there  Is  no  ''first  term"  which  can  he  nre-multioli ed.  Furthermore, 
the  index  of  the  bottom  quotient  used  must  be  orefixed  et  a  larve  enou»»h  value  to 
allow  for  the  worst  case  In  the  interval  allowed.*  Th< a  is  wasteful  <n  time,  since 
the  worst  case  will  seldom  be  encountered.  However,  for  a  function  like  tan"^x, 
which  is  so  troublesome  to  handle  with  series,  the  advantage  of  the  larger  Intervr.l 
allowed  by  a  continued  fraction  outweif^hs  all  its  disadvantages. 

Mecessary  Conditions  For  Efficient  Use  of  Series  In  the  604  or  605 

hhat  is  required  of  a  series  in  order  to  be  able  to  propram  its  com¬ 
plete  evaluation  within  a  604  or  605,  might  be  sunined  up  a  little  more  oreci  selv 
as  follows: 

There  must  exist  a  useful  interval  (a,b)  and  a  form  of  t^  such  that, 
for  X  in  (a,b), 

/  ’ 

(1)  j  ~  Kn,»)  is  recursive  and  is  both  sufficiently  ^mple 

and  apiall  enough  for  manipulation  within  the  comouter; 


*Soaetime8  aeries  arc  evaluated  this  way  -  i.e. ,  as  a  polynomial  of  high  dagree  - 
the  "starting*  (the  highest)  n  being  fed  in  with  the  arpument  end  subject  to  the 
discretion  of  the  c^er.  There  appears  to  be  little  advantage  to  this  met.hod.  tt  ■ 
takes  so  little  more  time  to  attain  maximum  accuracy,  once  the  series  Is  started, 
that  there  seems  to  be  no  reason  why  a  coder  should  have  to  compute  or  guess  at  the 
proper  n,  even  assuming  that  ha  knows  what  the  argument  wi 11  he. 

Those  familiar  with  CPC-II  timing,  will  realize  that  once  a  delay  cycle  is 
taken,  all  of  tt  is  then  available  for  coaoui^tion,  wh^ch  Is  several  times  the 
interval  allotad  to  a  card  cvcle  as  "compute  time." 
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(2)  an  N  can  be  chosen  which  will  keen  r(n,x)  within  necessarv 
bounds  for  n<  N  and  ».hich  w’ll  allow,  for  desired  d, 


(x)  4  10“P;  and 

N 


(3)  10*^  for  which  tp(x)  4  lOP  for  all  n,  5  for  wi  1  ch  ^  f: 


rnd  max  |^|  »  |^|  '  are  of  comparable  mapnitude  and  coToa.tible 
wi  th  p  in  ( 2) . 

It  is  best  not  to  be  too  stringent  vaith  (a,b)  or  some  usei’ul  series,  such  as  the 

lop  series,  are  likelv  to  be  summarilv  dismissed.  Likewise,  whether  or  not  r(n,x) 

is  "suf ficientlv  simnle"  deoends  nartlv  on  the  in^enuitv  of  the  nro^ra.mner.  The 

interdependence  of  the  demands  makes  more  preci .  statements  difficult. 

Examples  of  Fvaluatinr  Non-FI ementary  Functions  in  the  604. 

In  conclusion,  two  examples  are  riven  to  show  how  a  seH  •“s  can  sor-ie- 

ti.mes  be  used  for  rapid  evaluation,  within  restrictnri  intervals,  of  non-el  e.mnnt a rv 

functions  on  the  604:  aX 

(a)  the  probfibilitv  integral,  0(x)  =  dt,  and 

f7r/2 

(B)  the  elliptic  integral,  tCk)  =  /  (1-k^  sin^  0)2d0. 

J  0 

(a)  Since  0(x)  is  an  even  function  and  the  scries  for  e^  i  s  uniformly 
convergent  over  any  interval,  we  can  write,  where  k  =  '^j2lTT  , 

0(x)  =  k/  (1  -  +  ti_  _  +  .  .  .  +  (_l)n  .  .)  dt 

Jo  ?  2^21  p33.  pn„. 


=  k  X  - 

2-3 


'!-5  ?^-?!-7 


n  y2n+l 

+  .  .  .  +  (-1)  - 

2"-n'(2ntl) 


Lettinr  y  =  and  factoring  out  m  x,  we  have 

00 

0(x)  =  kx  X] 

o  (2n+l)n! 

For  X  in  [j3,b),  if  we  wish  to  carry  onlv  one  digit  for  QyJ  ,  v.e  must  have  ii  <  .'..472. 
To  keep  I  t^j<10  also,  b  is  further  restricted.  We  can  find  this  valic  bv  sottlnr 


kx 


2n+l 
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n  =  1,  2,  ...  and  solving  =  ^^2n+l)n\  ~ 

since  the  initial  kx  must  be  carried  through  all  the  terms  in  thp  actual  prograjnming 


to  save  storage.  Since  (x)  >  0,  for  all  n  and  x  >  0,  b  =  ml 


n  Jx 

L"J 


[■  i  s  obviously 


'n 


4  10  for  X  in  |^,b)  and  all  n.  Computati 


on 


necessary  and  sufficient  to  ensure 
gives: 

x^  =  4.2209,  xp  =  3.4676,  xj  =  3.2922,  X4  =  'i.71LB,  X5  =  3.313Q,  X(^  =”^.1798. 

Since  X  <  X-,  f.(x)  >  fc(x,),  i.e.,  r(4)  <  1  for  x  =  x  .  fince  the  series  is  an 
4  5  4  4  5  4  4 


exponential  type,  r(4+p)  ■4.  1  for  x  =  x^  and  hence  ^4(^4)  >  ^4+p(^4)>  P  ”  ^ 


f  •  •  •  • 


J'xl. 


Therefore,  x^  =  min.|^x^|.,  jf  we  similarly  compute  f  such  that 
in 


■'n 


4  100  for  X 


||b,b*),  we  find  that  b'  is  only  4.0326.  If  we  are  more  interested  in  accuracy, 
it  is  better  to  be  satisfied  with  b,  since  x  =  b’  is  a  pretty  largp  arnjment  for 
|^(x)  anytway. 

fte  can  carry  8  decimals  for  the  sum  by  ignoring  the  whole  numbers  in 
the  leading  terms,  for,  since  0  4  0(x)  <  1  for  nil  x,  these  must  all  cancel,  rnd 
the  worst  error  we  can  commit  by  ignoring  them  is  to  come  out  short  by  exactly  1. 

We  can  simply  add  1  to  the  result  and  ignore  whole  numbers  when  rendinr  out  the 
answer.  A  difficulty  arises  in  the  coraoutntion,  however,  since  r(n-l)  =  ^ 


(2n+l 


and  we  need  to  keep  9  decimals  for  (2^-l)t^  ^y  to  maintain  7  good  decimals  for  (^(x) 
Eut  if  n^  23,  2n^  +  n  >  1000,  which  is  too  larre  for  the  division  in  the  604 
without  using  extra  programs  which  are  not  available.  Limiting  n  to  22  reduces  b 


to  about  3.15  if  we  are  to  require 


'n 


to  become  less  than  10 


-7 


A  board  has 


been  wired  which  carries  7  decimals  for  t^  if  0£.  x  4.  3,  and  which  carries  6 
if  3£;  X  4  4.  For  4  4.  X  4  10,  the  computation  is  suppressed  and  a  constant 
.9999900  is  emitted  for  0(x);  also  a  negative  balance  selector  is  picked  up  to  be 

used  as  desired.  For  0  4  x  4.  3,  the  maximum  error  in  0(x)  is  lO”"^;  for  3£  x  4  4, 
it  is  10"^.  Speed  is  somewhat  better  than  33  ^cards  ner  minute  on  the  average. 
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B.  To  express  E(k)  as  a  series,  we  first  expand  (1-t^)^  in  a 
Taylor’s  series.  This  series  is  uniformly  convergent  for  t?  5-  b  ^  1;  although 
setting  t^  =  k^sin^0,  where  k  =  sin  0,  allows  t^  to  equal  1,  we  can  restrict  k^ 
to  jjDjb),  b  ^  1  to  bound  t^  away  from  1.  This  restriction  is  less  severe  than  the 
one  we  shall  eventually  have  to  make  anyway.  Vie  then  have, 

TT/p  TT/p  r-  TO 

E(k)  =£  (1-E2.in2^)id0  =  J'  j^l  - 


2 

kf^Mn^ 

7r\ 

Integrating  term  by  term  and  using  the  definite  integral 

7/2 


d(Z!. 


sin^n^  d0  =  1 . 3«  5« « ♦  ( JQ!  » 

2.4*6. ..{PnT  2 


TO  ^ 


we  have  E(k)=  ^jl-il£  -  ^ 


1 


00 


P,-  -  ^ r v>^--(.2n-i)r kfi.  LirZtn, 
2  I  4  ^  L  2.4.6...(2ny  J  2n-lJ  2  0 


where  t^  =  1,  t^^  =  -ik  ,  and  t^^^  =  "  2.^* 

This  is  not  too  difficult  to  program  on  the  604,  but  the  interval  is  restricted. 

It  will  probably  not  be  possible  to  allow  more  than  7  decimals  for  t  because  of 
the  magnitude  of  the  expressions  4n2-l  and  (n+l)2,  which  will  hav^  to  V'P  develooed 
as  a  multiplier  and  a  divisor,  respectively.  Seven  decimals  can  be  carried  for 
the  sum  which  allows  ^/7  to  be  inputted  at  read  time  and  carried  as  a  factor 
throughout  all  the  terms. 

For  k  =  .96296  =  sin  74°21’  +  ,11  4  lO""^,  but  this  is  not  Jin  uoper 

I  ICXJ 

bound  to  the  error,  since  the  series  is  non-alternating.  Fince  ^ 


4  as 


n 


00 ,  if 


=  k2/4  <  1/8,  we  can  say  the  error  is  no  greater  tiian  the  last 
term  included.  This  requires  k  <  .707107.  Actually  6  good  decimals  can 

probably  be  obtained  for  considerably  higher  values,  say  .8  or  even  .85.  In  anv 
event,  as  soon  as  is  less  then  the  allowed  tolerance,  the  comoutntion  will  ston 
if  it  has  been  programmed  properly. 


I 
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